Uncertainties of Antiproton and Positron Spectra from B/C Data and mSUGRA 

Contributions for Clumpy Halos 



A.M. Lionetto and A. Morselli 
INFN, Sezione di Roma II, via della Ricerca Scientifica,Roma, Italy and 
Dipartimento di Fisica, Universita di Roma "Tor Vergata", via della Ricerca Scientifica, Roma, Italy 



V. Zdravkovic 

Dipartimento di Fisica, Universita di Roma "Tor Vergata", via della Ricerca Scientifica, Roma, Italy 

We have studied the variation of e + and p top of the atmosphere spectra due to the parameters 
uncertainties of the Milky Way geometry, propagation models and cross sections. We have used 
the B/C data and Galprop code for propagation analysis. We have also derived the uncertainty 
bands for subFe/Fe ratio, H and He. Finally we have considered a neutralino induced component 
in the mSUGRA framework. We conclude that even in the case of diffusion convection model for 
the standard propagation, where secondary spectra gives a good fit of the data, SUSY contribution 
can not be excluded up to now due to the overall uncertainties. 
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I. PRODUCTION AND PROPAGATION OF COSMIC RAYS IN THE MILKY WAY 



o 
o 

CN| ■ We have chosen Galprop [l| as a public code for the treatment of propagation of all cosmic rays (CR) together. Our 
| scope has been to determine the total uncertainties in the calculation of e + and p top of the atmosphere spectra due 
to the uncertainties of geometrical and propagation parameters and cross sections. Here we give very short description 
| of processes included in propagation equation 
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where ip(r, p,t) is total phase space density. This equation is valid for all the types of particles. Isotropic diffusion 



is defined by the coefficient that depends from rigidity (momentum per unit of charge, p = ■§■) D xx = (3D (p/p ) 
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] inspired by Kolmogorov spectrum (S = 1/3) of weak magnetohydrodynamic turbulence In some models we have 
used a break in the index S at some rigidity po, with a value Si = below the reference rigidity pq. The convection 
velocity field V c , that corresponds to the Galactic wind, has a cylindrical symmetry and its z-component is the 
only one different from zero. It increases linearly with the distance z from the galactic plane, in agreement with 
magnetohydrodynamical models |3j. Reacceleration is determined by the diffusion coefficient for the impulse space 
D pp that is a function of the corresponding configuration space diffusion coefficient and of the Alfven velocity Va 
> ■ in the framework of quasi- linear MHD theory 0. Of course, Alfven velocity and convection velocity gradient in 
the Milky Way for reacceleration and convection terms are unknown parameters of propagation (there are no other 

• • ■ sources of information from which we could extract them, except the spectra of cosmic rays) and their possible range 
will be constrained by the analysis of fits of suitable data. The same procedure is valid for constraining the height of 
the Galactic halo and the other unknown parameters. This will be analyzed further in order to obtain all the possible 

$_i ' spectra of antiprotons and positrons using the sets of the constrained parameters. This procedure was already used 
in [13 for another propagation code, even for the antiprotons from neutralino annihilations. In the same reference 
can be found also the review of the propagation models in the Galaxy with an extended colection of historical and 
recent references. Injected spectra of all primary nuclei are power law in impulse dq(p)/dp oc p~ 7 . This power law 
approximation has been shown to be allowed in the framework of diffusive shock acceleration models, as well as a 
small break in the injection index 7 0, Q, @- Source term q(r,p) for secondaries contains cross sections for their 
production from progenitors on H and He targets. The last two terms in equationn&re loss terms with characteristic 
times for fragmentation and radioactive decay. Propagation equation is solved numerically, using the Crank-Nicholson 
algorithm. 

The heliospheric modulation in the vicinity of the Earth has to be taken in account. We have used a model in which 
transport equation (that describes diffusion processes in the hcliosphcrc and includes effects of heliospheric magnetic 
field and Solar wind) is solved in the force field approximation 9]. In this case Solar modulation is a function of just 
a single parameter that describes the strength of the modulation. All the dynamical processes are simulated with 
relatively simple changing of the interstellar spectra during the propagation inside the heliosphere, described by the 
formula 
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where U and p are energies and impulses of the interstellar and top of the atmosphere fluxes and cj> is the unique 
parameter that determines the solar modulation. 



II. UNCERTAINTIES OF CR SPECTRA 
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FIG. 1: Enveloping curves of all the good fits of B/C data for DR and DC model with their best fits inside and the best fit for 
DRB model. For the complete list of data see [13 ]. 



We were treating the two extreme cases of propagation models: the first that uses diffusion and reacceleration (DR) 
and the second that contains diffusion and convection (DC) [jj. Many parameters in the propagation equation are 
free and must be constrained by experimental data. Secondary to primary CR ratios are the most sensitive quantities 
on variation of the propagation parameters. The most accurately measured parameter is boron to carbon ratio (B/C). 
We have used a standard % 2 test: 



E 



, B/C e 

(er„ 



_($>B/C exp _ $ B/C tc ^2 
n ) 



(4) 



For DR model we have required reduced \ 2 l ess than 2 for the fit of the experimental data (fig. . We take 
the data with relatively small solar modulation between 325 MV and 600 MV where the force field approximation 
is better justified than for high modulation parameters. The allowed ranges of all parameters are given in table [I] 
Using them we have found the enveloping curves of all the e + and p spectra that present upper and lower bounds of 
uncertainty band. For e + relative uncertainty is varying from 30% under 1 GeV to 15% around 10 GeV (fig. El left) 
while for p is from about 10% up to 15%. For better visibility all the black and white figures that contain positron and 
antiproton data will contain just the part of the experimental data chosen to cover maximally all the energy range, 
while all the data are presented on the color figures. All data are taken from references [Pol llfij . 

For DC model (fig. 01 we have taken all the reduced % 2 values less than 2.8 for the variation of Do, diffusion indexes 
S\, below, and 5%, above the reference rigidity po = 4 GV, z, V c and injection index for primary nuclei 71 below the 



TABLE I: Allowed values for DR model propagation parameters. 



par./ Veil. 


z [/epe] 


Do [cm 2 s _1 ] 


s 




V A ^ ] 


minimal 


3.0 


5.2 10 28 


0.25 


2.35 


22 


best fit 


4.0 


5.8 10 28 


0.29 


2.47 


26 


maximal 


5.0 


6.7 10 28 


0.36 


2.52 


35 



TABLE II: Allowed values for the propagation parameters for DC model. 



par./val. 


z[kpc] 


Do[^f] 


82 


dVQ r km 1 
dz L skpc J 


7i 


72 


minimal 


3.0 


2.3 10 28 


0.48 


5.0 


2.42 


2.14 


best fit 


4.0 


2.5 10 28 


0.55 


6.0 


2.48 


2.20 


maximal 


5.0 


2.7 10 28 


0.62 


7.0 


2.50 


2.22 



reference rigidity pg = 20 GV and 72 above it. Positive variations around <5i = gave unsatisfactory fit. In order to 
take the smallest possible break of this index we have decided not to take negative 61 values. Allowed values for the 
propagation parameters can be found in table [H] The same analysis as for DR model gives relative e + error between 
20% above the maximum and 30% below it (fig. El left) while for p is about 20% around 20 MeV, 17% around the 
maximum and 25% around 20 GeV (fig. El right). We calculated PAMELA expectations for e+ (fig. EJ) and p (fig. El 
(parameters of the best B/C fit) using its geometrical factor and detector characteristics during the three years 
mission in which it will measure with high statistics various cosmic rays spectra. We have found also spectra that 
correspond to the parameters of the best fit of B/C data for subFe/Fe ratio (another important ratio for testing the 
parameters of the propagation models, fig. OJ, protons, He and e~ as well as corresponding uncertainties. For DC 
model fits are good, while DR overestimates p, He and e~. 




FIG. 2: Upper and lower bounds of positron (left) and antiproton (right) spectra due to uncertainties of the pro pagation 
parameters for DC, DR DRB model and DC and DRB model respectively. Experimental data are taken from (lia. Ha|. 

We have considered also the DR model with a break in the injection index for the primary nuclei spectra taken at 
rigidity of 10 GV 0,0]- We determined allowed values of the propagation parameters (tabic Hll)l demanding the same 
X 2 as for DR model (fig.^l. The positron and antiproton uncertainties are presented in fig. El Even if e + are fitted a 
little bit better at low energies and also primaries are fitted better, all of them remain overestimated while p spectra 
remain practically unchanged. For this reason on fig. Eland fig. El on the right, we refer just the DRB case. 

We have also seen how the obtained antiproton spectra change on variation of the most important antiproton 
production cross sections. Those are reactions that include all the types of hydrogen and helium. Antiprotons 
are created in the interactions of primary cosmic rays (protons and other nuclei) of sufficiently high energies with 
interstellar gas. Dominant processes are interactions of high energy primary protons with hydrogen, for example 
p + p — > p + p + p + p. Parameterization of this cross sections used in our version of Galprop code is give n by Tan and 
Ng ll2j. Other used cross sections, those of primary protons with other nuclei, are studied in reference [13j . From these, 
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FIG. 3: Left: Ratios (Sc+Ti+V)/Fe that correspond to the parameters of propagation for DC, DR and DRB model that gave 
the best fits of B/C data. The upper and the lower limits of all the graphics obtained with the parameters that give good fits 
of boron to carbon ratio are given for DC and DR model. Experimental data are taken from|l8|. Right: Top of the atmosphere 
spectra of e~ that correspond to the parameters of the best B/C fit for DC, DR and DRB model (lower curves) and local 
interstellar spectra (upper curves). Experimental data are taken from (Ifij . 

TABLE III: Allowed values for the propagation parameters for DR model with the break. 



par./val. 


z[kpc] 


D [^f] 


5 


7i 


72 




minimal 


3.5 


5.9 10 28 


0.28 


1.88 


2.36 


25 


best fit 


4.0 


6.1 10 28 


0.34 


1.92 


2.42 


32 


maximal 


4.5 


6.3 10 28 


0.36 


2.02 


2.50 


33 



the most important are those that involve helium and they contribute less than 20% of the total production of all the 
antiprotons. All the heavier nuclei together give just a few percents of the total production. Simultaneous settings 
of all the production cross sections to the maximum/minimum rise/lower the upper/lower parameter uncertainty 
bounds. Errors obtained in this way give contributions to the total uncertainties varying from 20% up to 25% in 
the case of DR model and, almost the same, from 20% up to 24% for DC model, depending of the energy range of 
the spectra. Non production cross sections, the so called tertiary component, correspond to inelastically scattered 
secondaries p + X — > p + X. Those processes tend to bring down the energies of the antiprotons of relatively high 
energies, flattening like that the spectra. But, even if their uncertainty is relatively big, this does not give relevant 
change of the spectra because tertiary contribution is very small. 

Changing of another exactly unknown quantity, He/H ratio, in a reasonable range from 0.08 to 0.11 gives a relatively 
small contribution, that vary from 3% - 7% depending of the energy, for both e + and p uncertainty, and for both of 
the models. Total uncertainties of e + and p are presented at fig. 0] and fig. |S] respectively. They vary from 35% up to 
55% for antiprotons and from 20% up to 40% for positrons roughly for both of the models in the current experimental 
data energy range. 

III. ANTIPROTON FLUX FROM NEUTRALINO ANNIHILATIONS 

In this section we take into account the possibility of a neutralino induced component in the p flux. Our analysis 
is performed in the well known mSUGRA framework [19] with the usual gaugino mass universality at the grand 
unification scale Mqjjt- The five input parameters arc defined as follows: 

m 1 / 2 , mo, sign(n), Aq and tan/3 , 

where mo is the common scalar mass, m 1 f 2 is the common gaugino mass and Aq is the proportionality factor between 
the supersymmetry breaking trilincar couplings and the Yukawa couplings, tan/3 denotes the ratio of the vacuum 
expectation values of the two neutral components of the SU(2) Higgs doublet, while the Higgs mixing fi is determined 
(up to a sign) by imposing the Electro-Weak Symmetry Breaking (EWSB) conditions at the weak scale. The pa- 
rameters at the weak energy scale are determined by the evolution of those at the unification scale, according to the 
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FIG. 4: Total uncertainties of e + fluxes and spectra that correspond to the parameters of the best B/C fit for DC (solid 
lines around the best fit) and DRB model (dashed lines around the best fit). Experimental data (from 16]) vs. PAMELA 
expectations for DC model. 

renormalization group equations (RGEs) 20]. For this purpose, we have made use of the ISASUGRA RGE package 
in the ISAJET 7.64 software [2l]]. After fixing the five mSUGRA parameters at the unification scale, we extract from 
the ISASUGRA output the weak-scale supcrsymmetric mass spectrum and the relative mixings. Cases in which the 
lightest neutralino is not the lightest supersymmetric particle or there is no radiative EWSB are disregarded. We have 
assumed a small clump scenario for the dark matter halo in our galaxy. By hypothesis the clump is a spherical 
symmetric compact object with mass M c i and density profile p c i (r c i)- We denote with / the dark matter fraction 
concentrated in clumps and we introduce the dimensionless parameter d 

d= 1 f d 3 r cl [ Pcl {r cl )} 2 
Po J <Pr c \p c \ (f c i) 

that gives the overdensity due to a clump with respect to the local halo density po = p(tq), where r$ is our distance 
from the Galactic center. In a smooth halo scenario the total neutralino induced p flux calculated for r = r (l is given 
by 



$ p (r ,T) = (a aan v)^^Bf (j^j C plop (T) 



(6) 



where T is the p kinetic energy, a ann v is the total annihilation cross section times the relative velocity, m x is the 
neutralino mass, B' and dN-> /dT, respectively, the branching ratio and the number of p produced in each annihilation 
channel / per unit energy and C prop (T) is a function entirely determined by the propagation model. 
In the presence of many small clumps the p flux is given by 

<S>f umpy (r ,T) = fd-<S> p (r ,T) (7) 

For the smooth profile we have assumed a Navarro, Frenck and White profile (NFW) [24| • 

The primary contribution to the p flux has been computed using the public code DarkSUSY 25] . We have modified 
the p propagation in order to be consistent with the DC model as implemented in Galprop. We assumed diffusion 
coefficient spectra used in Galprop with our best fit values for the diffusion constants Dq and S. In DarkSUSY the 
convection velocity field is constant in the upper and lower Galactic hemispheres (with opposite signs, and so it suffers 
unnatural discontinuity in the Galactic plane) while Galprop uses magnetohydrodynamically induced model in which 
component of velocity field along the Galactic latitude (the only one different from zero) increases linearly with the 
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FIG. 5: Total uncertainty of p fluxes and spectra that correspond to the parameters of the best B/C fit for DC (solid lines) 
and DRB model (dashed lines). Experimental data (from [l5|l vs. PAMELA expectations for DC and DRB model. 

galactic latitude 0. We have assumed an averaged convection velocity calculated from the Galactic plane up to the 
Galactic halo height z. 

We present here four different susy contribution to the p flux for different neutralino masses (obtained from a 
particular choice of mSUGRA parameters) and different dumpiness factor fd. Higher neutralino masses improve 
high energy data fit but only with the increase of the dumpiness factor because of the dependence from the inverse 
neutralino mass squared m x in the p flux formula We remark that for those models we have not computed the 
neutralino relic density. 

IV. CONCLUSIONS 

For e + in DR model, even when the uncertainties are included, the curve of the minimal e + production still remains 
above the experimental results. Breaking the primary spectra gently improves just the low energy part of the spectra 
below the maximum, but even the minimal predictions still remain bigger than the experimental results around the 
maximum as well as below it. On the other side, this slightly changes the best B/C fit (fig. P>. This is also not 
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FIG. 6: Different susy contribution to the total p flux with the DC model background. The dashed lines correspond to the 
background contribution, the punctuated lines correspond to the neutralino induced contribution while the solid lines correspond 
to total contribution 



sufficient to match protons and helium, that are still overestimated. In the end, e~ remain largely overproduced at 
low energies, even more than without the break. Uncertainty bands of p in DR models touch the experimental data 
from below. In this case the fit can be improved easily with any primary component that is coming from eventual 
neutralino annihilation or some other exotic contribution. For DC model all the results are in good agreement with 
the data, with only some problems with B/C data. 

As a possibility for the future considerations we would like to emphasize that it is natural to expect models that 
include both of the processes - convection and reacceleration, and this perhaps can be a key to solve simultaneously 
problems with positrons, electrons and B/C. But, we have not yet found such a kind of models. In any case, further 
measurements of p and e + spectra, primary to secondary CR ratios and Solar effects on them, as well as the precise 
determination of important nuclear cross sections seem to be crucial for answering the question of the models quality. 

In the framework of DC model exotic contributions remain possible at high energies (E > 20 GeV), and also they are 
not excluded at lower energies due to the relatively large uncertainties (if the Standard Model productions are smaller 
than the best fits). In the case of DR model background it is possible to improve the p data fit adding a neutralino 
induced component for a large part of the parameter space. Finally we have shown how a possible neutralino induced 



component in a clumpy halo scenario could give a contribution to the total p flux, in the case of DC model background. 
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